Transport on a Lattice with Dynamical Defects 
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Many transport processes in nature take place on substrates, often considered as unidimensional 
lanes. These unidimensional substrates are typically non-static: affected by a fluctuating environ- 
ment, they can undergo conformational changes. This is particularly true in biological cells, where 
the state of the substrate is often coupled to the active motion of macromolecular complexes, such 
as motor proteins on microtubules or ribosomes on mRNAs, causing new interesting phenomena. 
Inspired by biological processes such as protein synthesis by ribosomes and motor protein transport, 
we introduce the concept of localized dynamical sites coupled to a driven lattice gas dynamics. We 
investigate the phenomenology of transport in the presence of dynamical defects and find a novel 
regime characterized by an intermittent current and subject to severe finite-size effects. Our results 
demonstrate the impact of the regulatory role of the dynamical defects in transport, not only in 
biology but also in more general contexts. 



Spatial organisation and order is a fundamental char- 
acteristic exhibited by many complex systems. The high 
degree of coordination and organisation, particularly in 
living systems, is often achieved by forms of molecular 
transport, also crucial for material biosynthesis [T]. 

Transport in several biological processes is carried out 
on complex substrates of relative small sizes, such as mi- 
crotubule filaments or mRNA strands pp. These sub- 
strates are not rigid, and can either actively or passively 
change their conformation in time. Remarkably, the dy- 
namical features of the substrates are mostly coupled and 
regulated by their interplay with the transport process. 
In the case of protein synthesis, for example, ribosomes 
can encounter regions of the mRNA strand that obstruct 
the path and interfere with their movement. Those local 
regions act like a dynamical switch that obstruct, when 
folded, or allow, when unfolded, the passage of ribosomes. 
Importantly their folding-unfolding dynamics is strongly 
coupled to the presence of ribosomes: a secondary struc- 
ture cannot fold until the ribosome has moved forward. 

The traffic dynamics is highly enriched due to the cou- 
pling between the state of the substrate and the parti- 
cles, and despite its primary importance to the physics 
of transport phenomena, it has been usually neglected 
in the literature, with the exception of some particular 
cases related to microtubules dynamics [2]. Importantly, 
the phenomenology goes beyond of biological contexts; 
substrate dynamics and their bottlenecks are indeed cou- 
pled with the flow of carriers in many other disparate 
systems, such as intelligent traffic lights [3] or pedestrian 
traffic [1] . An extensive theory capable of studying the 
emergent effects of coupling between substrate and par- 
ticles transport is therefore necessary. 

In this Letter we present and study transport on a lat- 



tice with dynamical defects, i.e. sites whose features are 
dynamically coupled with the presence of particles. To 
do that, we use an approach based on a paradigmatic 
model in non-equilibrium statistical physics, namely the 
totally asymmetric simple exclusion process (TASEP) [5] . 
We define as dynamical defects sites with two possible 
conformational states that depend on their particle oc- 
cupancy. In the presence of dynamical defects, a new 
phenomenology appears. The current-density relation 
presents a plateau characterised by an intermittent flow 
of particles that cannot be estimated by standard mean- 
field arguments. Remarkably, in this new regime we find 
large finite-size effects that strongly affect the transport 
characteristics. 

The Model - The standard TASEP is a model of parti- 
cles moving unidirectionally on a one-dimensional dis- 
crete lattice with a fixed hopping rate 7. Steric in- 
teraction excludes that more than one particle can oc- 
cupy the same site, and overtaking is not allowed. The 
average TASEP current of particles J as a function 
of particle density p is given by the parabolic relation 
Jtasep(p) = Tp(l — p) in the large lattice limit. The 
presence of one or more defects in the lattice, defined as 
"slow" sites with hopping rates smaller than 7, leads to 
a reduction of the maximal current flow [5] and confers 
a truncated parabola shape to the current-density rela- 
tion J(p). The defects considered so far in literature, 
however, are static since hopping rates associated to the 
slow sites are constant. 

The periodic boundary case illustrates well the main 
phenomenology of the model. We consider a ring of L 
sites on which N < L particles are allowed to hop in one 
direction, fixing the density at p = N/L. A region of size 
d represents the dynamical defect and has an intrinsic 
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FIG. 1. (a) Mechanism for a 1-site dynamical defect: the 
site s (grey square) closes or opens with rates / and u, re- 
spectively. Particles hopping at rate 7 cannot enter closed 
regions. Moreover, if a particle occupies the site s, the clos- 
ing of that site is forbidden, (b) Sketch of the phase diagram 
of the system in the f/<y and 11/7 space. The different regimes 
are separated by dashed lines. 



two-state dynamics coupled to the presence of particles. 
The defect can pass from the closed state to the open 
state with rate u. The inverse process occurs with rate 
/ only when all d sites of the defect are empty. Our 
results indicate that the single-site defect d—1, Fig.[lji, 
exhibits the main features of the model. Therefore in this 
Letter we focus, for sake of simplicity, on this case. The 
model proposed here can be applied to describe further 
systems; for instance, junctions in transport networks [7] 
can be thought of as particular dynamical defects with 
influx dependent dynamics. 

Numerical Simulations - We have performed 
continuous-time Monte Carlo simulations based on the 
Gillespie algorithm |S]. We have studied the model in a 
large parameter space taking several lattice sizes L (up 
to 4000) and varying the rates 7, u, / in order to explore 
all the different dynamical phases sketched in Fig. [T]b. 
Similarly to other cases involving static defects or lo- 
calized blockages in ID systems, our model shows a 
reduction of the current with respect to the homoge- 
neous TASEP; the current-density relation is a trun- 
cated parabola with a constant plateau value, see Fig. 
SI 0. In large systems, when the opening rate u de- 
creases with respect to the other typical time scales, the 
current-density plateau lowers to smaller and smaller val- 
ues, occupying a wider interval of densities and merging 
with the TASEP parabola 7/d(l — p) only at the extremi- 
ties (for p — > and p — > 1). As in [6 , clusters of particles 
form before the defect so that the average density profiles 
are usually characterized by a sharp separation between 
a High Density (HD) and a Low Density (LD) phase. 
Depending on the relationship among the parameters 7, u 
and / we obtain different dynamical regimes, for which 
we determine the characteristic timescales via the compu- 
tation of the probability distribution function of the time 



lag between the passage of two consecutive particles at a 
given site [9]. Each regime exhibits a different slope in a 
log-linear scale defining characteristic timescales. Three 
main dynamical regimes appear: (i) TASEP-like regime: 
when u > f > 7 or u > 7 > /, due to the rapid open- 
ing, the current is not affected by the defect. The system 
therefore behaves like a homogeneous TASEP and the 
only relevant timescale is 1 /j (white region in Fig. [IJd) ; 
(ii) Static defect-like regime: when / > u > 7 the dy- 
namical defect acts like a static one, allowing the passage 
of particles in a way that can be described with a single 
effective rate q < 7 (top-right gray region in Fig.JTJ)), sim- 
ilar to the model in [6]; (iii) Intermittent regime: when 
7 > / > u there are two strongly separated timescales, in 
contrast to the former cases. The short timescale 1/7 is 
the time separation between the passage of two consec- 
utive particles during an opening event, while the long 
timescale 1/u corresponds to the time intervals during 
which the defect is closed and obstructs the passage of 
particles (black region in Fig. [if}). This is a novel regime 
caused by the presence of the dynamical defect. The re- 
maining cases (7 > u > f and / > 7 > u) are crossovers 
regimes between the intermittent and the homogeneous 
or static-defect TASEP-like behaviors, respectively. 
Strikingly, numerical simulations also show that the sys- 
tem is subject to severe finite-size effects. If not otherwise 
specified, the following results are intended to be valid in 
the limit of large systems. We will come back to the 
finite-size effects in the last section. 

Mean-field approaches - To compute the current- 
density relationship we use a finite-segment mean-field 
(FSMF) approach [TO]: the dynamics within the defect 
is treated exactly while the injection and exit rates in 
and from the dynamical defect are given by effective val- 
ues. By imposing the current continuity and assuming 
that the system is split in HD and LD phases separated 
by a dynamical defect at site s, the densities have to be 
coupled as phd — 1 — Pld = 1 — Ps {ps denotes the den- 
sity of the defect site and we assume that the site s is in 
the LD phase) . We then model both HD and LD parts as 
open TASEPs with effective entry and exit rates [7] , with 
a 1 = 7(1 — p s ) being the injection rate and /3' = 7(1 — p s ) 
the depletion rate to and from the defect, respectively. 

We solve the related master equation in the steady- 
state [9J and compute, as a function of p s ,f,u and 7, 
the probability to find a particle on site s, which is by 
definition equal to the normalized density p s . One then 
obtains an expression for p s as a function of all other 
parameters [9i . The plateau current is then given by 

J "plateau 

= 7p s (l - p s )- This procedure can be extended 
to larger defects (d > 1). 

To validate the FSMF theory, in Fig. [2^, we show the 
relative difference A J/ J = (J — Jfsmf)/J between sim- 
ulations and the FSMF in all the different regimes of 
the phase diagram. Data are taken in the middle of 
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FIG. 2. (color online) (a) Relative difference AJ/J — (J — 
Jfsmf) I J between simulations and FSMF approach. Cir- 
cles and diamonds show the crossover from the intermittent 
regime to the TASEP-like regime, while squares and trian- 
gles show the crossover from an intermediate-non intermittent 
regime to the slow-site like regime, see Fig. lb. (b) Average 
closing times (tf) (in seconds) for different densities. Vertical 
dashed lines represent the boundaries between the different 
regimes crossed (see, dotted-dashed line in the inset). 



the plateau of the current-density relation (p = 0.5), 
also to avoid deviations due to finite-size effects, which 
will be discussed in the last section. This analysis re- 
veals that the FSMF approximation is not appropriate in 
the intermittent regime (circles and diamonds in Fig. [2^,, 
u/f < 1), while it gives reasonably good results in the 
other regimes, in particular in the TASEP-like and the 
slow-site-like regimes. 

Intermittent Dynamics- We now turn our attention 
to the intermittent regime, where the FSMF fails; we 
start by analysing the average time between consecutive 
opening and closing events (tf), Fig. ^p. When the cou- 
pling between the conformation of the defect and the 
presence of particles is weak, i.e. u is the largest rate, 
the average time scales as (tf) ~ f~ x . But when the 
coupling is strong, i.e. 7 is the largest rate, the former 
expression is corrected as (tf) = [/(l — /5)] -1 , where p 
is the probability to find a particle on the defect site 
given that the site s is open. This probability can be 
approximated in some limiting cases: in the TASEP-like 
regime, u is the fastest rate and hence the defect is al- 
most always open. In this case p ~ p, and the current 
is very close to the one predicted by the pure TASEP. 
In contrast, in the intermittent regime J(p) exhibits a 
plateau [5 , within which the opening-closing dynamics 
is independent of the total density. In this case particles 
are blocked for long times behind the closed defect that 
occasionally opens, then allowing for a collective passage 
of particles, as confirmed by kymographs or space-time 
snapshots of the system [S|. Just after the opening of 
the defect, p can be estimated from a simple mean-field 
approach: dp/dt = 7p/(l — p) — 7/5(1 — p r ), where pi and 
p r denote the density of particles before and after the 
defect, respectively. Approximating p\ by 1 and p r by 0, 
one gets p = 0.5 in accordance with numerical simula- 



tions (Fig. [2)d). 

Hence, an intermittent dynamics is established so that 
no net flow passes through the defect for a time of order 
\ju and then a large current J open — 7/5(1 — p) flows for 
a time (tf). This gives an average plateau current 



JiMF = 7/5(1 - p) 



u/f 



1-p + u/f 



(1) 



which we refer to as intermittent mean-field current 
(iMF). Figure [3^, shows that deep in the intermittent 
regime (where p ~ 0.5) this approach provides a substan- 
tial improvement to predict the current in the plateau, 
compared to the FSMF. 

Finite-size effects- As anticipated, the system presents 
strong finite-size effects, remarkably pronounced in the 
intermittent regime; here the current-density relation of 
small lattices becomes asymmetric (Fig. [3^,) showing a 
reduction of the current for small densities and enhance- 
ment for large densities. Figure [3^, shows that the value 
of J in the plateau is not affected by the system size. 
However, the plateau disappears when systems with just 
a few sites are considered (Fig. PJd) . 



J 



0.5 



x L = 200 
O 500 


1000 
- FSMF 


iMF 


ac^g^ c^c 


00/ 








X 

1 


1 1 


1 



J 



o°°o 

or 



- Ye(i-e) 

iMF 

O u = 0.07 „ 
x 0.05 O/ x 
+ 0.02 <?■•■ 
• 0.01 9- 

+ 



<y ¥ 

7V 



H 



0.2 0.4 0.( 

Q 



0.8 1.0 
a) 



0.5 

Q 



1.0 
b) 



FIG. 3. (color online) (a) In the intermittent regime truncated 
parabolas are found only in the large L limit: smaller systems 
show current reduction (enhancement) at low (high) densities. 
The FSMF calculation is unable to predict the plateau value, 
while the iMF does. The parameter used are / = ls _1 ,u = 
0.01s -1 and 7 = lOOs^ 1 . (b) For very small systems (here 
L — 16) the iMF (dashed lines) correctly approximates the 
numerical J(p) relation (/ = Is -1 , 7 = 100s _1 ). 



Such a behavior can be rationalized as follows: while 
in the thermodynamic limit the probability to find a par- 
ticle on the defect site, when open, is well approximated 
by p ~ 0.5, in very small systems (e.g. L = 16 in Fig.^) 
the unstable HD-LD interface moving from the dynam- 
ical defect quickly relaxes over the whole system to the 
homogeneous TASEP density. Therefore, when the de- 
fect allows the passage of particles, all the sites can be 
considered to be identical and the density on the defect 
p ~ p. Thus, identifying p with p in Eq. (nl) gives the 
quantitative good agreement shown in Fig. [3p. More- 
over, this reasoning allows determining the lower bound- 
ary for the Jl(p) relation for very small system sizes L, 
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the upper boundary Joo (p) being the truncated parabola 
profile predicted by the iMF theory. For intermediate 
system sizes L the situation is again different. Remark- 
ably, the current for p > 0.5 is enhanced compared to the 
the current obtained in the thermodynamic limit, while 
it is reduced for p < 0.5 (Fig. [3^). This effect can be 
understood by analyzing the relaxation dynamics after 
an opening event by integrating the system of L coupled 
ODEs describing the occupancy of each site in the lat- 
tice 0. 

Intermittence and related finite-size effects have strong 
consequences on the density profiles too. While, similar 
to the static-defect case, outside the intermittent regime 
there is a sharp phase separation between the HD-LD 
before and after the dynamical defect (Fig. [4^,), the pres- 
ence of intermittence induces relevant boundary effects in 
the density profiles (Fig. [4)d) . This is a signature of the 
transient relaxation of the HD-LD interface during the 
time (tf). Such strong correlations are associated to the 
non-stationary dynamics of finite systems, and are par- 
ticularly evident at densities for which current reduction 
(or enhancement) occurs. 
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FIG. 4. (color online) Density profiles for a lattice with L = 
500,/ = Is" 1 , it = 0.01s -1 . (a) For 7 = Is' 1 there is a clear 
HD/LD separation with a sharp domain wall, in agreement 
with the FSMF (horizontal dotted lines), (b) In presence of 
intermittence (7 = 100s _1 ) longer correlations are established 
and simulations deviate from the FSMF prediction. 

Conclusion - In this Letter we have introduced the 
concept of transport on a lattice in the presence of lo- 
cal interactions between particles and the substrate, also 
referred to as dynamical defects. This concept is key to 
understand natural transport phenomena occurring on 
substrates with a fluctuating environment. It allows un- 
derstanding fundamental biological processes such as pro- 
tein synthesis or intracellular traffic, where the local con- 
formation of the substrate can deeply influence the char- 
acteristics of the flow of molecular motors. The obtained 
results, however, are general and therefore applicable to 
other transport processes, such as vehicular or human 
traffic, or synthetic molecular devices. 

The phenomenology is presented and studied by 
means of the totally asymmetric simple exclusion pro- 



cess (TASEP), a prototypic model of transport in non- 
equilibrium physics. In this framework we have provided 
original mean-field arguments that allow reproducing and 
rationalizing the rich phenomenology of the model. We 
have thus unveiled a novel dynamical regime, character- 
ized by an intermittent current of particles and induced 
by the local interaction and competition between par- 
ticle motion and the defect dynamics. Importantly, we 
have found that a particle-lattice interaction triggers se- 
vere finite-size effects that have a counterintuitive strong 
impact on transport. For different system densities, the 
small size of the lattice unexpectedly reduces or enhances 
the flow of particles, inducing an asymmetric current- 
density profile. 

Given the small size of biological substrates, the 
physics of the new intermittent regime is highly rele- 
vant to understand protein synthesis and motor protein 
transport. Bursts of gene expression have been often re- 
ported and our results provide a physical mechanism 
to explain the contribution of the translation process to 
them. Moreover, intermittent behavior can strongly in- 
fluence motor protein current fluctuations that can be 
measured in state of the art experiments |12j . In the 
biological context of mRNA translation, our results on 
finite-size effects would correspond to an increase in pro- 
tein production for small mRNA strands compared to 
longer strands; interestingly, highly expressed proteins 
constituents of ribosomes are short [T5]. 

The model presented here can be extended by in- 
cluding, for example, larger interaction sites or differ- 
ent boundary conditions. However, the observed phe- 
nomenology does not change: if the defect is extended 
(d > 1) the finite-size effects are even more pronounced 
than in the d = 1 case [14] , and lattices with open bound- 
aries, more common than lattices with periodic bound- 
ary conditions in practical applications, present the same 
characteristics, although edge effects can be relevant if 
the defect is moved close to the boundaries [T4"] . 
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